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Abstract 

If the vacuum manifold of a field theory has the appropriate topological structure, the theory 
admits topological structures analogous to the D-branes of string theory, in which defects of one 
dimension terminate on other defects of higher dimension. The shapes of such defects are analyzed 

o ' 

CN ■ numerically, with special attention paid to the intersection regions. Walls (co-dimension 1 branes) 

terminating on other walls, global strings (co-dimension 2 branes) and local strings (including gauge 



fields) terminating on walls are all considered. Connections to supersymmetric field theories, string 
theory and condensed matter systems are pointed out. 
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I. INTRODUCTION. 



Topological defects arise whenever spontaneous symmetry breaking leads to a topologi- 
cally nontrivial manifold of vacua or ground states. Defect configurations involving spatially 
nonuniform order fields are of considerable interest in any field theory which possesses them 
as their topological quantum numbers distinguish them from vacua with uniform order fields. 
A complete understanding of the physical significance of defects requires a good knowledge 
of their field and energy profile and much effort has been devoted over the years to this 
subject (see □D for reviews). Applications range all the way from fundamental physics to 
cosmology and condensed matter phase transitions. 

For a symmetry breaking phase transition G — > H, the vacuum manifold is given by 
the coset space M = G / H , and p-dimensional defects are classified by the homotopy group 
7r d-p-i(Jvt)> where d is the spatial dimensionality. For field theories with a relatively simple 
field content and abelian homotopy groups the associated topological defects are essentially 
isolated fundamental objects - they do not intersect. 

In more complicated field theories new possibilities arise involving composite defects or 
intersecting arrays of defects. Such configurations add a new layer of mathematical and 
physical richness. The simplest examples are o-dimensional defects whose boundaries are 



themselves (p — l)-dimensional defects 



s 



Another example is furnished by theories 



with non-abelian homotopy groups, such as the quaternionic fundamental group of RP 2 , in 
which topological defects of a fixed dimension (strings) intersect. Some time ago Carroll and 







Trodden 13, |8J considered the case of defects which terminate when they intersect defects 
of equal or higher dimensionality, such as strings (p — 1) ending on domain walls (p = 2). 
The defects on which other defects end were termed Dirichlet topological defects (DTD), in 



analogy with D-branes in string theory 



ii 
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141 ]. extended objects on which 



fundamental strings can terminate. It is important, of course, to keep in mind that there 



are important differences between the two sets of objects, as emphasized in [21, |lf 

Related configurations are of interest in several contexts. Strings terminating on walls 



can arise in Yang-Mills theories 



18 



y 



20 



21 ] as well as in grand unified models, where 
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171 ] . In nonabelian gauge theories, 



they may be of use in tackling the monopole problem 
non-intercommuting cosmic strings provide a further example. In addition, in higher spatial 
dimensions, the required symmetry breaking schemes can lead to a nonabelian unbroken 
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gauge group, and the phenomenon of confinement for the resulting defects. The rich variety 
of symmetry breakings and associated ground state manifolds in condensed matter systems 
provides a particularly concrete context to study topological defects, including the Dirichlet 
class studied here. In fact an N = 1 vortex filament in the nonchiral superfluid 3 He — B 
phase can terminate at a domain wall interface separating this phase from a superfluid 
3 He — A phase [2], thus realizing a string terminating on a wall. 

Since interesting nonlinearities arise both in the cores of such defects and at the points 
at which distinct classes of defects intersect, it is important to have explicit solutions at 
hand to study, rather than simply a classification of the associated vacuum structure which 
guarantees their existence. Some work along these lines has already been performed, par- 



ticularly in the case of defects networks in supersymmetric theories 
related work in the context of D-branes see 



34, 



35 



4J |2 



29 



30 



3lj (for 



36|). This is the focus of the 



current paper, whose outline is now provided. In Section |H] we treat the case of walls ending 
on walls, followed in Section ITTT1 bv the case of global strings ending on walls. In Section HVl 
we add gauge fields and treat the case of local strings terminating on walls. In each case 
we give explicit field and energy profiles. Throughout this paper we work in three spatial 
dimensions. 



II. WALLS ENDING ON WALLS 



The Model 



The simplest example of a DTD is a wall ending on another wall. Walls are produced by 
the spontaneous breaking of discrete symmetries. Consider a theory with three fields, 0, ipi 
and invariant under three Z 2 symmetries. Such configurations are a special case of Z v 
walls 



•po, invariant under three Z 2 symmetries. Such configurations are a special case of Z„ 
[37], and are related to solutions in supersymmetric and supergravity jillQl 



theories. 
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with Lagrangian density 



7] 



L = \d ll ct> &><j> + \ d^ x d^i + \ d^ 2 d*i> 2 - V(<p, </>i, ih) , (2) 
where we choose the most general renormalizable potential invariant under ZP x zf x Zf 



V(<t>, ^ 2 ) = X^ 2 - v 2 ) 2 + X^j M + ^ - w 2 + g{<p 2 - v 2 )} 2 + - WW - V 2 ) • (3) 

For = 0, the minima of V are the 8 points given by 



±v and < 



(^i,^ 2 ) = (±*,0) 
or 

(V^Va) = (0,=hS) 



(4) 



Since the potential is positive-definite in this case, the presence of h requires either ipi or ip 2 
to vanish, after which the first term in the potential is minimized by <p = ±t> and the second 
term by ip\ = ±w or ip 2 = ±w. For fi > the degeneracy of the 8 vacua is broken, leaving 
only 4 minima of the potential, at 



(V'ljV'a) 



(±u;,0) 
(O,±i0) 



(5) 
(6) 



where v is the solution of the cubic equation 



IX^X^v +QX^g/j,v - (8X<pXipv + // ) v — 2X^/i (gv +w) = 



(7) 



satisfying lim v = v , and w is given by 



w 



w 2 + g(v 2 - v 2 ) + 



2X,i, 



1/2 



(8) 



The equations of motion for static fields are then 

V 2 = 4A 0(0 2 - v 2 ) + AgX^l + $ - * 2 + s(0 2 - v 2 )} - //(^ - ife) 

VVi = 4A^ift^ + ipl -w 2 + g(<f) 2 - v 2 )} + 2h^l - 

V 2 ^2 = 4A^ 2 [^ 2 + - * 2 + ^ - *> 2 )] + 2^ 2 i/>2 + 2/^0^2 ■ 



(9) 
(10) 

(11) 
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This model is constructed to admit solutions consisting of domain walls terminating 
on other domain walls. Such a configuration is intrinsically two-dimensional, and hence 
is independent of one spatial coordinate, which we choose to be z. What remains is a 
boundary value problem for an elliptic system of partial differential equations, for which we 
must prescribe the value of the fields on the boundary. 

B. The Numerical Problem - Boundary Conditions 

Our choice of boundary conditions is illustrated in Fig. Q so that we have Dirichlet 

© © 




x 
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FIG. 1: Boundary conditions for walls on walls. The circles (o) represent Dirichlet boundary 
conditions and solid lines represents Neumann boundary conditions. 

boundary conditions 






%y = 


—00) 


= —V 
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%y = 


—00) 


= , 


^2 


%y = 


— 00) 


= w , 



at y — > —00, and Neumann boundary conditions everywhere else except at the upper corners, 
at which we choose 





[x = 


±00. 


y = 


+00) 


= v , 




[x = 


±00, 


y = 


+00) 


= ±w , 


^2 


[x = 


±00. 


y = 


+00) 


= . 
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These choices define a wall in the ifj\ field, terminating on a wall in the field - the Dirichlet 
wall. 

The numerical algorithms we use require homogeneous boundary conditions and so we 
introduce new fields 

4> = Uq + F , 

V>i = u x + Fi , (14) 
V> 2 = u 2 + F 2 , 

and define 

n / \ r / \ l + tanh(2y) r , N , . 
■ro(ac, 2/) = 7os(y) H ^ L/ouW - Jo^U/ = +oo)J , 

Fi(x,y) = fu(y) tanh(2x) + 1 + t^hgy) _ ^ = +oo) tanh(2a;)] ? (15) 

F 2 (x,y) = f 2s (y) ■ 

The functions fi U (x) and fi. s (y) are the solutions of the systems of ordinary differential 
equations obtained by restricting to the "up" boundary (y = +oo) and the "side" 

boundary (x = +oo) respectively, with i = 0,1,2 corresponding to 0, i/ji,i[>2 respectively. 
These functions are easily found using standard algorithms with boundary conditions 



(16) 





fou 


(±oo) 


=v , 






flu 


(±oo) 


= ± w , 






jW 


(±«)) 


=0 . 




/os( 


-oo) = 


— V 


/ Os (+oo) 


= I) 


fu( 


-oo) = 





/i s (+oo) 


= If 


f2s( 


-oo) = 


w 


/2 S ( + 00) 


= 



(17) 



The symmetries of our problem are such that the solutions on the boundary at x = — oo 
are then ip\(x = — oo, y) = —fu(y) and ip2{x = —oo,y) = / 2s (2/)- Note that one solution 
f2u(%) = is trivial. 

We are now in a position to solve numerically our partial differential equations in the 
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interior. We use a multigrid algorithm 43] to solve for the fields Ui, satisfying 

Ui = at the y = — oo boundary and at (x, y) = (±00, +00) , 
d x Ui = at x = ±00 , (18) 
dyUi = at y — +00 . 

C. Results - Walls on Dirichlet Walls 

We present our results as a set of plots of field profiles and energy densities, providing a 
number of ways of viewing configuration shapes. We have already defined the field profiles, 
and the energy density for static fields is given by 

The field domain wall configuration is plotted in Fig. El The intersection with the ipi 
wall at x = is clearly visible. The corresponding configuration for is shown in Fig. |3J 




FIG. 2: The <fi Dirichlet wall configuration. 

Note the domain wall for y = +00. Finally the ip2 field configuration is plotted in Fig. 0] 
It has the form of a wall with half the height of the Dirichlet wall in the sense that it goes 
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FIG. 3: Field configuration for 




FIG. 4: Walls on walls field configuration for ip2- 

from (at y — +oo) to w — 1 (at y — — oo). The energy density plotted in Fig. illustrates 
the merging of the walls at the origin. A corresponding contour plot for the energy density 
in Fig. El shows the wall merging very clearly. 
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FIG. 5: Energy density for walls ending on walls. 
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FIG. 6: Contour plot of the energy density for walls ending on walls. 



III. GLOBAL STRINGS ENDING ON WALLS 



A. The Model 



In the simplest case, strings are generated by a breakdown of a U(l) symmetry. In 
analogy with our approach in the previous section, consider the following transformations 



9 



on one real scalar field and two complex scalar fields and ip 2 . 

Z 2 ■ {<P ~<f>, ipi ^2) 

: ^ _, e --i (20) 

C/(l) 2 : ^ 2 ^ e"^ 2 ^2 • 

In this section, both U(l) symmetries are global, i.e. u>i is space-time independent. The 
Lagrangian density is 

L = \ <9 M + <9^i c^i + 9^2 «9^2 - V(<j>, ihM , (21) 

with potential 

v(<j ) ^ 1 ^ 2 ) = x^<f ) 2 -v 2 ) 2 + o^r+i^r-^+s^ 2 -^)] 2 

+ %ii 2 i^r-Mi^r-i^i 2 ) • (22) 

Consider the following static ansatz 

(f) = 0(r, z) 

i/> 1 = R 1 (r,z)e ie (23) 
^2 = "02 (r, 2;) , 

where 0, i?i and ^2 are real functions, for which the equations of motion are 

V 2 = 4A 0(0 2 - v 2 ) + AX^Rl + 4, 2 2 -w 2 + g^ 2 - v 2 )} - ^{R\ - ^) , (24) 
V 2 Rt = 2\^ J R 1 [R\ + ip 2 -w 2 + g(<p 2 - v 2 )} + hR 1 ip 2 - y,<\)R x , (25) 
V 2 ^ 2 = 2A^2 [R\ + 4% - w 2 + g(<j) 2 - v 2 )} + hR\i\) 2 + ^2 . (26) 

B. The Numerical Problem - Boundary Conditions 

We now proceed as in the previous section, imposing the boundary conditions illustrated 
in Fig. 

We choose Dirichlet boundary conditions for the bottom boundary 



01 




Z = 


—00 ) 


= -v , 




Rv 


( r > 


z = 


—00) 


= 0, 


(27) 


^2' 


( r > 


z = 


—00) 


= w , 
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FIG. 7: Boundary conditions for a string ending on a wall. The circles (o) represent Dirichlet 
boundary conditions, the squares (□) represent mixed boundary conditions (Dirichlet for some 
fields, Neumann for the remaining ones) and solid lines represent Neumann boundary conditions. 

and at the upper corner we impose 



0(r = 


oo, 


z 


= +oo) 


= V 


Ri(r = 


oo, 


z 


= +oo) 


= w 


ip 2 (r = 


oo, 


z 


= +oo) 


= 



(28) 



The boundary r = requires a more detailed discussion. Taylor expanding the fields 
around r = 0, we have 



= ^2a n (z)r n , R 1 = ^b n {z) 

n n 

Equations [2~Hl26j . in the limit r — > 0, become 



^2 = ^C n {z) 



" ' ' ' + 4a 2 (z) + = fo[a (z), b (z), c (z)} 



r dz 2 
b (z) d 2 b (z) 



3b 2 (z 
cAz) 



+ 



+ 4:C 2 (Z) + 



dz 2 
d 2 c (z) 



r dz 2 
Demanding regularity at the origin gives 



f 1 [a (z),b (z),c (z)] , 
f 2 [a (z),b (z),c (z)} . 



(29) 



(30) 



ai(z) 


= 


dr<f>\ r=0 


= o , 


bo(z) 


= 


R l\r=Q 


= , 


ci(z) 


= 


d M r=0 


= . 



(31) 
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For the other boundaries we imposed Neumann boundary conditions for all the fields. 

As for the case of walls on walls, we define new fields with homogeneous boundary con- 
ditions. Defining the "up" boundary to be at z = +00 and the "side" boundary to be 
at r = +00 we again introduce the functions fi U {r) and fi S (z) which are solutions to the 
appropriate systems of ordinary differential equations with boundary conditions 

d r fou(0) = / Ou (+oo) = v , 

A«(0) = /i„(+oo) = w , (32) 
d r f2u(0) = / 2u (+oo) = , 

/Oa(-Oo) = -V / Os (+Oo) = V , 

/la(-Oo) = /la(+Oo) = W , (33) 

/2a(-Oo) = W /2s( + 00) = . 

Again we have a trivial solution /2 U (r) = 0. Following the previous section, we define the 
following functions 

F (r, z) = f 0s (z) + 1 + ta 2 nh(2z) [f 0u (r) - f 0s (z = +00)] , 

F ± (r, z) = f ls (z) tanh(r) + 1 + ^h(2z) _ f ^ = ^ ^ ^ 

F 2 (r,z)=f 2s (z) , 
and introduce new fields 

cf) =u + F , 

^ = Ul + F 1 , (35) 
ip 2 = u 2 + F 2 , 

obeying homogeneous boundary conditions 

Ui — at the z = —00 boundary, at (r, z) = (+00, +00) 
and for u\ at the core of the string, 
d r ui = at the r = +00 boundary and for i — 0, 2 at the core of the string, 
d z Ui = at the z = +00 boundary . 
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C. Results - Global Strings Ending on Dirichlet Walls 



To help in understanding our plots, note that the energy density for static fields in this 
section is given by 

2 -i / n i\ 2 / n / \ 2 / n ; \ 2 



*-* .—K£H(S + (£W£j + (36) 

The field domain wall configuration is plotted in Fig. |H1 The intersection with the ipi 
string at r = is clearly visible. The corresponding configuration for ipi is shown in Fig. |§1 




FIG. 8: Field <j), the Dirichlet wall. 

Note the string for z = +00. Finally the ip2 field configuration is plotted in Fig. It has 
the form of a wall (except for an interaction with the string for r = and z > 0) with half 
the height of the Dirichlet wall in the sense that it goes from (at z — +00) to w — 1 (at 
z = —00). 

The energy density plotted in Fig. ^2 illustrates the merging of the string with the wall at 
the origin. A corresponding contour plot for the energy density in Fig. ^] shows the string 
ending on the wall very clearly. 
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FIG. 9: Field ipi, the global string. 




FIG. 10: Field V>2 for a global string and a wall. 
IV. GAUGE STRINGS ON WALLS 

More interesting connections with physical models can be made by introducing gauge 
fields. In this section we take the first steps towards this by finding numerical solutions for 
local cosmic strings terminating on Dirichlet domain walls. 
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FIG. 11: Energy density for a global string ending on a wall. 




FIG. 12: Contour plot of the energy density for a global string ending on a wall. 
V. THE MODEL 

We modify the Lagrangian density of the last section by introducing a U(l) gauge field 
Ap, under which the field ipi is charged. 

£ = \ d,<p 8^ + D^i D^ipi + d^ 2 d^ 2 - \ F^F^ - V(<t>, M 2 ) , (37) 
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where 

DfM = <9 M — ieA^ . (38) 

The equations of motion for and ip2 are unchanged from the case of the global string, since 
we have constructed the model so that the relevant winding occurs for the ipi field only. 
Thus, the modified equation of motion for ip\ is 

dV 

d^ifr - ied^A^x) + _ - ieA„{dyi - teA^) = . (39) 
Imposing the Lorentz condition d^A 1 * = for a static gauge field yields 

V • A = , (40) 

and equation (|3U|) then becomes 

8V -> ^ ^ 

VVi = 77— -2ieA- V^ + e 2 {A- A)^ . (41) 
dipt 

We are interested in the ansatz 

<p = <f>(r,z) , 

ip 1 =R 1 (r,z)e ie , (42) 
^2 = ifa(r,z) , 

where 0, i?i and ip2 are real functions. As is well known for cosmic strings, the imaginary 
part of equation (jUj) gives 

^+^ = 0, (43) 

with solution 

A = A6, (44) 

so that the Lorentz condition becomes dgA = 0. The equation for R\ then becomes 
„ „ d 2 R l 1 <9i? R x d 2 R 1 



AR 

2A v ,i?i [i^ + - * 2 + #(> 2 - £ 2 )] + kRi^i - A*0-Ri + 2e + e 2 A 2 R 1 , (45) 



and the equation for the gauge field is 

UA a = f = 2eIm[^i(<9 Q - ieA )^} , (46) 

which, for our ansatz, is 

(V 2 A) e = 2eR 2 1 (r~ 1 + eA) . (47) 
16 



A. The Numerical Problem - Boundary Conditions 



The geometry of the problem is identical to that for global strings, and so we refer the 
reader again to Fig. [7| We choose Dirichlet boundary conditions for the bottom boundary 



0(r, z = - 


-oo) = 




Ri{r,z = - 


-oo) = 


0, 


^2(r,y = - 


-oo) = 


w , 


A(r,y = - 


-oo) = 


0, 


b(r = oo, z 


= +00 


)=v 



(48) 



and at the upper corner 

Rx(t — 00, z — +00) = w , 
ip 2 ( r = 00, z = +00) = , 
A{r = 00, z = +00) = . 

At the core of the string we choose 

d r (j)(r — 0, z) = d r ip 2 (r — 0, z) = 
R x ( r = 0,z) = A(r = 0,z) = , 



(49) 



(50) 



and for the boundary at r = +00 we impose Neumann conditions for all the fields. The 
choice of these boundary conditions defines a string ending on a Dirichlet wall. 

We again solve ordinary differential equations on the boundaries at z — +00 and r = +00. 
We keep the same notation as in the previous section, but now include the gauge field A by 
defining f^ u {r) and f^ a {z) with boundary conditions 

hu(0) = <9 r / 3 u(+oo) = , (51) 



/ 3s (-oo)=0 / 3s (+oo) = 0. (52) 
Once again f2 U ( r ) = and f^z) = are trivial solutions. 
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As before we define 

F (r, z) = f 0s (z) + 1 + ta 2 nh(2z) [fou(r) - f 0s (z = +00)] , 

F 1 (r, z) = f u (z) tanh(r) + 1 + tanh ^) [/lu(r) _ fu{z = +oo) tanh(r)] 

2 - (53) 

F 2 (r,z) = f 2s (z) , 

„ , . 1 + tanh(2,2) „ . , 
Fa(r,z) = ^-LJ. faa (r) , 

and new fields 

4> = u + F , 
^ = u x + Fx , 

(54) 

i/j 2 = u 2 + F 2 , 
A = m 3 + F 3 , 

with homogeneous boundary conditions 

Ui = at the z = —00 boundary, at the corner (r, z) = (+00, +00) 
and, for i — 1, 3, at the core of the string, 
<9 r Wj = at the r = +00 boundary and, for i — 0, 2, at the core of the string, 
<9 2 Wj = at the z = +00 boundary . 

B. Results - Local Strings Ending on Dirichlet Walls 

We plot our results as before, with the energy density now given by 



(55) 

dr J \r J \ dz J 2 

We also plot the magnetic field associated with the gauge field of the string, given by 



1 + Rl (- + eA\ + ( ^ V + l -B 2 + V(<f>, R lt fa) 



- cM A 1 d(rA) A 

5 = -_r + -^-^z . (56) 

The field domain wall configuration is plotted in Fig. The intersection with the ipi 
string at r = is clearly visible. The corresponding configuration for ip\ is shown in Fig. 
Note the string for z = +00. The field ip 2 configuration is plotted in Fig.E3 It has the form 

18 



FIG. 13: The <fi field configuration: Dirichlet wall for local strings. 




FIG. 14: The local string field 

of a wall (except for an interaction with the string for r = and z > 0) with half the height 
of the Dirichlet wall in the sense that it goes from (at z — +oo) to w — 1 (at z — — oo). 
Finally the gauge field configuration is shown in Fig. We can verify that A — > r~ l as 
(r, z) — > +oo. We plotted the r component of the magnetic field, B r , in Fig. [T7| The z 
component of the magnetic field, B z , is illustrated in Fig. ^] A 2-D plot of the vector field 
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FIG. 15: Field configuration ip2 for a local string ending on a wall. 




FIG. 16: The gauge field A. 

B is shown in Fig. EH The energy density plotted in Fig. 120 illustrates the string merging 
at the origin with the Dirichlet wall. A corresponding contour plot of the energy density in 
Fig. |^ shows the string ending on the wall very clearly. 
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FIG. 17: Component B r of the magnetic field. 




FIG. 18: Component B z of the magnetic field. 
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FIG. 19: 2D plot of the magnetic field B. 




FIG. 20: Energy density for a local string ending on a wall. 
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FIG. 21: Contour plot for the energy density for a local string ending on a wall. 
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